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1 . Introduction 

Processor sharing (PS) is a mathematically tractable approxi- 
mation to time sharing, a procedure followed in manv actual 
computer systems. In effect, PS assigns to each gob of one i, 

( i = 1»2,...) present :or processing l/’ith of th e total Drcces s ing 
efrort; ecuivalently , a single icb with Markovian service rate 
u completes processing m (c,r-dt) with probability (g./i)dt + o(dt). 
One advantage or PS is tnat short jobs are not tr acred behind 
long gobs, as is possible in a FC-FS discipline. 

Various mathematical results have been obtained aoout certain 
processor sharing models. An early example was the caper of 
Coffman e t al . (19 70) . Recently extensive results have been ob- 
tained for Markovian systems by D. Mitra (1931), and for non- 
Markovian single-server Poisson arrival systems by T. Ctt (1934) 
and V. Ramaswami (1934). 

This paper is a continuation of work reported in Gaver, Jacobs, 
and Latouche (1984), henceforth GJL, where emphasis was placed on 
proposing and evaluating simple approximations to the distribution 
of delay experienced by a particular "tagged" gob approaching a 
time-shared processor. In that paper it was shown that under heavy 
traffic conditions, or if the tagged job duration became large, 
then the distribution of tagged job response time (also called 
sojourn time by others) approaches the normal or Gaussian 
distribution . 
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We first show that the analysis of our system and the problem 
solution can naturally and conveniently be conducted in work time 
rather than in ordinary clock time. Subsequently we focus upon 
heavy-traffic approximations to the distribution of response time, 
R (T) , given the actual work (computer time) requirement, T. 
Specific models proposed and examined are (i) for a system of 
many identical terminals that independently submit jobs or pro- 
grams according to the same Markovian process, and (ii) for a 
system having two terminal types, each of which submits jobs in 
a manner governed by its own Markov process. The methodology 
extends to general k terminal types as well, and to other models. 

The approximation solutions are evaluated for accuracy by 
means of Monte Carlo simulations. 
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2 . 



The Work-Time Concept 



Imagine that a tagged job requiring T units of processing 
approaches the computer. Assume that it arrives when the system 
is in steady state. After that initial moment it undergoes 
processing, at various rates governed by the amount of its 
accompaniment, until T units of service or work are accumulated, 
at which point it departs after a random delay of R(T) . In 
following the tagged job's delay, it turns out to be convenient to 
measure time in terms of the amount of actual work or processing 
that has been accomplished on the tagged job . Thus let {X(w) , 
w >_ 0 } denote the number of programs or jobs undergoing service 
at a moment when exactly w units of processing have been accomplished 
on the tagged job. The instantaneous rate of accrual of clock or 
response time at work time w is clearly X(w) : if X(w) = 1 then 

the tagged job is alone and response (clock) time and work time 
advance at the same rate, while if X(w) = 17 the tagged job is 
accompanied by 16 others and 17 units of response time accrue for 
every single work time unit. It follows that 

T 

R ( T) = / X ( w) dw . (2.1) 

0 

For the models to be considered, the process X(w) is a birth- 
and-death, or simple Markov, process related to N(t), the number 
of jobs in the system at clock time t; its transition rates, 
after adjusting for tagged job entry, are seen to be 

A^dt = A (N-i) dt = A (N-i)idw = A^dw 

( 2 . 2 ) 

t w 

y^dt = ydt = y i dw = y^dw 
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for 1 <_ i £ N where N is the total number of terminals in the 
system. The term idw is required to allow the work-time process 
to advance appropriately in clock time. Henceforth we drop 
the superscripts, allowing the context to imply the appropriate 
rate . 

The approach taken here invokes the work time concept 
described above to facilitate calculations, first for a single 
terminal or job type situation, but later on for a system with 
more diversified traffic patterns. 
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3 . Heavy Traffic Analysis of a Single Terminal Type Markovian 

System 

Suppose that N terminals have access to a single computer. 

Each terminal has Markovian demand rate A , and expected service 
time y ^ . Service times are assumed to be independently and 
exponentially distributed. The discipline at the computer is 
PS. It is clear that if N(t) is the number of terminals that 
have submitted jobs that are undergoing service at the computer 
at (clock) time t, then (N(t) , t >_ 0 } is a Markov process in 
continuous time that is identical to the classical single repairman 
problem; see Feller (1957) , p. 416. This is so, since if 
N(t) = i > 0, then each individual job or program receives (dt/i) 
units of processing time in (t,t+dt) , and hence departs with 
probability y(dt/i) + o(dt) , but the probability that some job 
departs is iy(dt/i) + o(dt) = ydt + o(dt) . It has been shown by 
Iglehart (1965) and by Burman (1979) that under heavy traffic con- 
ditions (N -*■ «) one may approximate N(t) by a suitable Gaussian 
process, namely the Ornstein-Uhlenbeck process. This fact alone 
enables one to study the distribution of R(T) , and to deduce 
approximate normality; see GJL for a first analysis. 

3 . 1 Diffusion Approximation in Work Time . 

Here is a diffusion process approximation for X(w) . On the 
basis of intuition write down the stochastic differential equation 
for X(w) , the approximation to X(w) : 

dX ( w) = A [N-X ( w) ] X (w) dw - yX(w)dw 

(3.1) 

+ /A [N-X (w) ] X (w) + yX(w) dB(w) 
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where (B(w), w >_ 0} is standard Brownian motion. The first 
right-hand-side term represents infinitesimal drift of X(w) , 
while the second is the diffusion or infinitesimal variance term, 
the form of which is obtained from the observation that arrival 
and departure processes compete like independent Poisson processes 
in short time periods. Now suppose that, as N °° , 



where m(w) is a deterministic function of time, and { Z (w) , 
w >_ 0 } is a stochastic process, the properties of which must be 
discovered. Substitute (3.2) into (3.1) to obtain 

Ndm(w) + /NdZ (w) = A [N-Nm(w) -/NZ (w) ] [Nm(w) + /NZ (w) ] dw 



+ /A [N-Nm(w) -/NZ (w) ] [Nm(w) + /NZ (w) ]+y [Nm(w) + /NZ (w) ] dB (w) . 

Next isolate terms of order N and /N; the result is, after 
stipulating that A' = AN, a constant as N + «=, 



X ( w) - Nm(w) + /N Z (,w) 



(3.2) 



- y [Nm(w) +/NZ (w) ] dw 



(3.3) 




A ' [ 1-m (w) ] m( w) - ym(w) , 



(3.4) 



0(/N) : dZ ( w) 



{A 1 [l-2m(w) ] - y } Z ( w ) dw 



+ /A ' [l-m(w) ]m(w) +ym(w) dB (w) ; 



(3.5) 



the stochastic differential equation (3.5) is of Ornstein- 



Uhlenbeck (O-U) form; see Arnold (1974). 
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Next obtain the approximate long-run mean as the solution 
of (3.4) with dm/dw = 0, examining only the heavy-traf f ic 
situation in which A ' > y : 

m(«) = 1 - t^t , A ' > y (3.6) 




If the above solution is used to define the stochastic 
differential equation parameters there results 

dZ(w) = ( - A ' +y ) Z (w) dw + /TJTd"- Ty7'A" r ) ) dB(w) , (3.7) 

which suggests that {Z(w)} can be considered an O-U process with 
constant coefficients; namely 

dZ (w) = -pZ(w)dw + odB(w) , (3.8) 

the solution to which is 

Z(w) = Z ( 0 ) e -pw +o / e~ pu dB ( u) . (3.9) 

0 

2 

The parameters p = (A'- y) and a = 2y(l-y/A'). 

3 . 2 Response Time Evaluation 
Let 

T T 

R(T) = / X(w)dw - / [Nm(w) + /NZ(w)]dw (3.10) 

0 0 
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approximate the response time; in this approximation R(T) is 
normally distributed (Gaussian). First, 



T 

E[R(T)] - E [R (T) ] - N / m(w)dw = N(l- 

0 




(3.11) 



Second, 



T 

Var[R(T)] - Var[R(T)] = var [/N j Z (w) dw] 

0 



(3.12) 



N{ E [ / Z (w) dw f Z(u)du)]-(E[ f Z ( w) dw ] ) ^ } ; 



0 



0 



0 



for ease of writing we have left the initial condition Z(0) 
implicit. In order to evaluate the above, recall that the tagged 
job approaches the server when the latter is in equilibrium, 
i.e., at t = 00 . It may be shown that the diffusion approximation 
for N ( t ) , the number undergoing service at clock time t, is 



where a(t) is a deterministic function of time and {Y (t) } is a 
particular Ornstein-Uhlenbeck process. A similar analysis to 
that leading to (3.6) and (3.7) yields 



N ( t) = Na ( t ) + /N Y ( t ) 



(3.13) 




(3.14) 




and 



Q 



E [ Y (oo ) ] = o, Var [Y (°° ) ] 



= JL = JL . 
NX " X' ' 



(3.15) 



see GJL for a derivation; (3.13) provides the initial condition 

for evaluating moments of R(T) , using (3.10) . Identify Z(0) , 

the initial value of the work time noise process Z (w) , with Y( ot ) . 

T 

According to (3.9) , this implies that E[ f Z(w)dw] =0. In order 

0 

to compute the Var [R(T) ] it is next necessary to evaluate the 
following integral: 



T T 

I (T) = E [ / Z (w) dw j Z(u)du] 

0 0 



T T 

= E[2 J Z(w)dw j E(Z(.u)|Z(w)]du] 

0 w 



r i , _ , 

= 2E[ / Z(w)dw j Z (w) e K U W dui 

0 



w 



(3.i6: 



T 

= 2E [ f Z (w) dwZ (w)i(l - e" p (T ' W) ) ] 

J n 



= 2E[ / E [Z (w) 2 | Z (0) ]-(l - e' p(T W) )dw] 



Square (3.9) and take the expectation to see that, condition- 
ally on Z ( 0 ) , 
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I(T) 



= 2E [ 



0 



/ (Z (0) 2 e 2pw + 0 



-2pw. .1., -p (T-w) . , 
2 ^(l-e M ) ) — ( 1-e M ) dw 



— E 1 Z ( 0 ) 2 J [ / e- 2ow (l-e-° <T - w) )dw] 
P o 



♦ 4 / T (l-e‘ 2pu ) (l-e- p<T - w) )dw 

P 0 



Now put E{2(0) 2 ] = E[Y(®) 2 ] = u/X 1 = a 2 /2p to see that 



I (T) = 4 / T H-e‘ p(T * w) Jdw - 4lT -id -e' oT ) 1 

P 2 0 P 2 0 



^[l-i(l-e- pT )] 



Thus it follows that 



Var [ R (T ) ] - Var [R( t) ] 



NT 2^11 --^(1 - e" pT ) ] 

P 



(3.14) 



To terms of order T this agrees with (4.10) of GJL; not surprisingly 
the additional factor in (3.14) can actually provide numerical 
results superior to those of GJL. 

The form of the heavy traffic approximation, namely the limit- 
ing normal form with parameters (3.11) and (3.14) , can be more 
rigorously validated by use of the theory of convergence of 
suitably normalized sequences of semigroups of transformations; 
see Burman (1979) . Details appear in an appendix to this paper. 



10 



Heavy Traffic Analysis of a K-Terminal-Type Processor 
Sharing System. 



4 . 



Consider the following natural extension of the previous 
model. The processor is jointly utilized by K sets of terminals, 
each generating distinctive job types. There are terminals 



respectively. Again the discipline at the computer is PS. Of 
course this is not the same as a situation in which all terminals 



terminal. The latter model can, however, be studied in an 
analogous heavy-traf f ic manner, as can other interesting models. 

4.1 A Diffusion Model for the Work-Time Process. 

Let (X^(w) , i = 1 , . . . , K} represent the number of jobs of all 
types present at the computer at work time w. The present model 
implies that {X^(w)} is a multivariate or vector-state birth and 
death Markov process. We choose to study a diffusion approxima- 
tion (X^(w)} to {X i (w) } that is described by the following system 
of s . d . e . : 



in the ith set, and arrival rate and service rate are A. and u. 
— i 




dX i (w) 



K ~ 

A . ( N . -X . ( w) ) ( J X (w) ) dw - u . X . (w) dw 

111 i t K 11 



111 




i ~ K - 

A . (N . -X. (w) ) y X (w) + ]i . X. (w) dB. (w) 
111 ^ 11 1 



= 1 



(4.1) 



i = 1,2, ... ,K 



where {B^(w)} are mutually independent standard Brownian motion 



11 



or Wiener processes. The work time process is a transformation 
of the clock time process; in particular, the drift of the i th 
component of the clock time process {N^(t)} is seen to be 

N. (t) dt 

A i (N i -N i (t) )dt ' y i k ' (4.2) 

N (t) 
k =1 K 

which exhibits the processor-sharing effect in the term multiply- 
ing y^. Multiplication by the total in service, ^N^(t), converts 
to the work time transition rates, in analogy with (2.2). 

Now once again approximate by writing 



X i (w) = N i m i (w) + /N i Z i (w) , i = 1,2 ,...,k; 



:4.3) 



m. (w) and {Z. (w) } are to be determined, subject to the normaliza- 
1 K 1 

tion N = £ N. 00 but with 

k =1 



N ./N 

l 



l . , 0 < i . < 1 , 

l i — 



(4 .4a) 



and 



NX. - X!, 0 < X ! < 00 , X .’ > y . . (4.4b) 

i l-i—ii 

Conditions (4.3) and (4.4) are referred to as the heavy traffic 
normalization (HTN) . The result of isolating terms according to 
order in (4.1) is : 
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0 ( N) : 



(4.5) 



dm . (w) 
1 



K 



dZT~ - A!(l-m.(w))j; £ k m k (w) -l^m. 

K — 1 



(w) 



0 (/N) : dZ i (w) 



K K 

= A! {-( l in. (w) ) Z. (w)+/T7(l-m. (w) ) £ /ITz, (w)}dw 
k=l 1 1 1 k=l K K 



- y . Z . (w) dw 
l l 



(4.6) 



K 



+ y A| (1-iru (w) ) £ (w) +u^nr (w) dB^(w) 

k=l 



for i = 1,2, ...,K. Thus (nu (w) ; w >_ 0 ; i = 1,2,...,K) must 
be found by solving a system of ordinary first-order, but non- 
linear differential equations, while (4.6) shows that 
{ Z . (w) ; w >_ 0 , i = 1,2, . . . , K } is a mul ti variate Ornstein- 
Uhlenbeck process. 

4.2 A Diffusion Model for the Clock Time Process. 

In order to provide the initial conditions encountered by 
the tagged job, it is necessary to study the clock- time process 
N^(t); see (3.13). The corresponding approximation has s.d.e. 



N.(t) 

dN . (t) = A . (N . - N . (t) ) dt - y . -=*-= dt 

1 111 1 x\ 



l (t) 



k=l 



(4.7) 



/ - V t} 

+ V A i (N i" N i (t) )+U i~K — Z dB ± (t) i 



- 1,2 K 



l N v (t) 



k=l 
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Now invoke the HTN : 



N. (t) = N . a . (t) + /N . Y . (t) , 

i 11 11 



K 



and again N = £ N, , with 

k=l K 



N . /N 

l 



i ■ , 0 < l . <1 , 

l — l — ' 



but 



U i /N 



u ! , 0 < y , A . > y . 

i — i i i 



The result of isolating terms is 



0 ( N) 



da^ ( t) 
dt 



A i (1 -*i (t)) - - 



( t) 



I Vk (t) 



k=l 



0(/N) 



Y, (t) 



dY i ( t ) 



-A i Y.(t)dt-u;{ — 



l Vk (t) 



k=l 



K 









■}dt 



( l 4 k a k (t)) 



k=l 



a (t) 

+ W^U-a.UD+ij! — dB i (t) 



l Vk 111 



k=l 



(4.8) 



(4.9) 



(4.10) 



(4.11) 



(4.12) 



These equations closely resemble those describing the work time 
approximation; again the semigroup approach is applicable. 

If a long-run solution to the 0(N) term exists in work time, 
and consequently dm^/dw -* 0 as w -*■ °° , the result is the system 
of equations for nr (») 5 nr: 



K 



Pi U -m i ) l 



k=l 



*k m k - m i 



= 0 



(4.13) 



where = A|/y^ = NA^/y^. Now these same equations are satis- 
fied by a presumed long-run solution in clock time, i.e., if 
da^/dt ^ 0 in (4.11); for a^(°°) = a^: 

K 

p i ( 1 - a f) " a j/ l ^ k a k = 0 • (4.14) 

k=l 



Consequently the long-run solutions in work and clock time agree 
at the 0 ( N ) term level; this means that the long-run mean number 
present in both clock and work time agree: 



E[N i (t)J - Na i (oo) = Nm i («) ~ EtX^t)] . (4.15) 



Next substitute these long-run results in the s.d.e. to see 

that as t,w Y. and Z. are essentially the same process. 

K 1 1 

Put S = y 2,, a, to simplify writing. Then 
k=l k k 



dY i (t) 



_A i Y i (t) dt ~ 



Y i ( t ) 



dt + 



A i (l-a i ) /2, i 



K 

l /ITY k (t) 
k=l K K 



dt 



+ /2A i (l-a i ) dB i , 



(4.16) 
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or 




+ /2 A . ( 1-a . ) S — dB . ( t) 

1 1 /7T 1 



/s 



( 4 . 17 ) 



and a direct comparison with the corresponding equation in work 
time, ( 4.6 ) , shows that the long-run behaviors of the two 
processes (Z^(w) } and {Y^(t) } are identical except for a constant 
time- scale change ; for large w and t, 



4.3 Response Time 

We discuss the response time under these conditions: a 

tagged job approaches the processor when the latter has been 
operating for some time, so the long-run clock time distribution 
prevails; after arrival, the job remains present until the total 
work time accumulated on the job is T, the requested service 
time, giving 



{Z i (w) } and { ( ^) } 



( 4 . 18 ) 



have the same probability law; i.e., 
finite-dimensional distributions and 
limiting distribution. 
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T K K T 

R(T) = / 1 [ X . ( w ) dw ] - l j X . ( w ) dw 

0 i=l 1 i=l 0 1 



(4.19) 



K 



T 



K 



T 



= N 1 l . / m . ( w ) dw + /N £ /FT j Z . ( w ) dw 

i=l 1 0 1 i=l 1 0 1 



where it is understood that the initial condition for the 



distribution from the clock time process. In view of (4.18) , 
this is equivalent to removing the initial condition by the 
long-run distribution of the work time process itself. 

Since the long-run situation is being discussed it is first 
necessary to solve the steady-state version of (4.5) 



Then the solution provides parameters for the long-run version 
of (4.6); here written in matrix form 



Z^(w) integrand in (4.19) is given by the approximate stationary 



K 



0 




y • m . , i = 1,2, ... ,K. (4 .20) 

ii' 



dZ (w) 



A Z (w) dw + a dB 



(4.21) 



Now to find the variance of R (T) , append the row 




K 

7 v / T~ Z . (w) dw 
i=l 1 1 



(4.22) 
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to the former drift matrix A of (4.21) , and consider the system 

dZ*(w) = A*Z*(w)dw + a*dB* (4.23) 

the solution to which can be formally written out in terms of the 
appropriate fundamental matrix, and computed in terms of eigen- 
values and eigenvectors of the matrix A*. See e.g. Arnold 
(1974), Chapter 8 and Coddington and Levinson (1955) for details. 
A convenient way of formalizing the calculations is actually by 
using Laplace transforms. Unfortunately, no truly simple 
formulas result. Finally, the covariance matrix C(w) of the 
components of Z_*(w) satisfies the matrix differential equation 

dC(w) 

— = (A*)C(w) + C(w) (A*) ' + (a) (a) ' 

where 1 denotes transpose; the initial conditions are provided 
by the long-run distribution in clock time, or in view of (4.18) , 
of the work time process {Z} itself. it is the K+l— diagonal 
element of C(w) , evaluated at w = T and multiplied by N that 
provides the required approximate Var(R(T)]. 
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5 . Simulation Studies of the Accuracy of the Normal Approxima - 
tions to the Distribution of Response Time 

In this section we use simulation to study the numerical 
accuracy of normal approximations to the distribution of the 
response time . Two continuous time Markov chain models were 
simulated. In one there is a single terminal type; in the 
second, a two-terminal type system is examined. Two normal 
approximations were evaluated: one results from a central 

limit theorem, and the other results from applying the previously 
derived diffusion approximation to the Markov processes. 

5.1 A Single Terminal Type Markovian System. 

Let X(w) denote the number of other jobs undergoing service at 
a moment when exactly w units of processing has been accom- 
plished on the tagged job for the single terminal type 
Markovian model of Section 3. Since {X(w) ; w >_ 0 } is a Markov 
process and 

T 

R (T) = / ( X (w) +1) dw , 

0 

it follows that there are constants m(c) and o (c) such that 

R (T) - m( c) T 

o ( c ) /T 

converges in distribution to a standard normal distribution 
as T - °° (cf. Keilson (1979) , p. 121) . Call this a central 
limit theorem (CLT) for such a process. In this case 
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m(c) 



1 + l TT ( j ) j 



where tt is the stationary distribution of {X(w); w >_ 0 } and 

a formula for evaluating o(c) is given in Keilson (19 79) . The 

CLT normal approximation states that R(T) has a normal dis- 

2 

tribution with mean m(c)T and variance a(c) T. The CLT mean 
is the true mean for R(T) under steady state (of. GJL) . 

The CLT normal approximation should be increasingly accurate 
as T becomes large, despite values of other system parameters, 
including the number of terminals. 

The derivation of the heavy traffic (or diffusion) approxi- 
mation is detailed in Section 3. In summary, the HT approxi- 
mation is that R ( T) has a normal distribution with mean 



N (1 





1 



and variance 

NT ~2 ^ 1 " ^ F (1 ' e " PT) ] ' 

P 

where 



2 

a 



2U(1 




and 



p — N A + y . 

The mean and variance of the HT approximation are easier to 
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compute than those for the CLT . It is anticipated that the HT 
approximation should be increasingly accurate as N becomes large 
when heavy traffic conditions prevail, i.e. ^ < 1 . it is 
inapplicable under other circumstances. We have conducted simu- 
lations to assess these anticipations. All simulations were 
carried out on an IBM 3033 computer at the Naval Postgraduate 
School using the LLRANDOMII random number generating package 
(see Lewis and Uribe (1981)). 

Conditional response times given the number of jobs being 
processed at the time of arrival of the tagged job were simu- 
lated; the tagged job required T time units of processing. For 
each initial condition, 500 replications were done. Sample 
moments and relative frequencies were computed for each initial 
condition giving conditional response time sample moments, and 
selected response time relative frequencies, i.e., estimated 
probabilities of response times in selected ranges. Uncondi- 
tional sample moments and relative frequencies were then computed 
by multiplying each conditional moment or relative frequency 
by the appropriate stationary probability of there being j jobs 
present at the time of arrival of the tagged job and then 
summing over all possible j . The stationary probability is 
of the form kAr(j) where k is chosen so that the probabilities 
sum to 1 (cf. Kelly (1979)). A detailed description of the 
simulation program can be found in Pornsuriya [1984] . 

Simulated and approximating means and standard deviations 
for various values of N, A, and u appear in Tables 1-3. Some 
discussion of specific cases now follows. When N = 10, A = 25 
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and y = 100, it follows from (3.6) that approximately 
(10) (1 - j^-) = 6 jobs are being processed along with the 
tagged job; thus the traffic is moderate in this case. When 
N = 10, A = 15, y = 100, then on the average 3.3 are being 
processed with the tagged job; a rather light traffic case. When 
N = 25, X = 10, y= 100, then on the average 15 jobs are being 
processed; again a moderate traffic case. The HT mean is 
lower than the simulated mean for N = 10 . For N = 25 it 
agrees with the simulated mean. As mentioned before, the CLT mean 
equals the true mean. The CLT standard deviation approaches 
the simulation value as T becomes large, as anticipated. Also 
as anticipated, the HT standard deviation is closer to the 
simulated one for the moderate traffic cases than for the 
light traffic case. In order to assess the degree of normality 
of the distribution of R(T) , the a-quantiles for each approxi- 
mating normal distribution were computed. The relative 
frequency of being less than or equal to each a-quantile was 
computed using the simulated data. The results appear in 
Tables 4-6. For N = 25, A = 10, y = 100, the HT approximation 
appears to describe the data well for all values of T. For 
N = 10, the HT approximation does better than the CLT for the 
moderate traffic case of A = 25. For the light traffic case, 

N = 10, A = 15, both approximations do poorly for small T = 0.01, 

the mean work request time. The CLT does well for 

large T = 0.10, which is 10 times the average service time. 

Note, however, that all simulations have been carried out 
for modest system sizes, N. If N grows to say 50, or 100, the HT 
approximations can be expected to improve correspondingly; 
they are often not bad even at the level of N = 25. 



TABLE 1 



Simulated Mean and Standard Deviation for R(T) 
and Their Approximating Values 



N = 10, A - 15, p = 100 



TIME T 




Mean 


Std. Dev. 


0.01 


Simulation 


.0404 
( .0002) * 


.0172 




CLT 


.0404 


.0323 




HT 


.0333 


.0238 


0.025 


Simulation 


. 1005 
( .0005) 


.0375 




CLT 


. 1010 


.0510 




HT 


.0833 


.0535 


0 .05 


Simulation 


.2016 
( .0010) 


.0622 




CLT 


.2019 


.0722 




HT 


.1667 


.0919 


0.10 


Simulation 


.4029 
( .0015) 


.0940 




CLT 


.4039 


.1020 




HT 


. 3333 


.1462 


* The 


standard error of 


the estimate 


of the mean 
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TABLE 2 



Simulated Mean and Standard Deviation for R(T) 
and Their Approximating Values 

N = 10, X = 25, u = 100 



TIME T 




Mean 


Std. Dev 


0.01 


Simulation 


.0606 
( .0002) * 


.0503 




CLT 


.0605 


.0245 




HT 


.0600 


.0160 


0.025 


Simul at ion 


.1507 
( .0004) 


.0315 




CLT 


.1513 


.0288 




HT 


. 1500 


.0314 


0.05 


Simulation 


. 3021 
( .0008) 


.0497 




CLT 


. 3027 


.0548 




HT 


. 3000 


.0481 


0.10 


Simulation 


.6036 
( .0012) 


.0748 




CLT 


.6053 


.0776 




HT 


. 6000 


.0706 



★ 

Standard Error for the Estimate of the mean 
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TABLE 3 



Simulated Mean and Standard Deviation for R(T) 
and Their Approximating Values 



N = 25, A = 10, u = 100 



TIME T 




Mean 


Std. Dev. 


0.01 


Simulation 


.1498 
( .0002) * 


.0812 




CLT 


.1500 


.0390 




HT 


.1500 


.0254 


0 .025 


Simulation 


. 3759 
( .0006) 


.0503 




CLT 


.3750 


.0617 




HT 


. 3750 


.0497 


0.05 


Simulation 


. 7506 
( .0010) 


.0804 




CLT 


. 7500 


.0872 




HT 


. 7500 


.0760 


0.0625 


Simulation 


.9381 
( .0011) 


.0909 




CLT 


.9375 


.0975 




HT 


.9375 


.0863 


* 

The 


standard error 


of the estimate 


of the mean 
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TABLE 4 



Simulated Probability (Relative Frequency) that the 
Response Time is Less than or Equal to the Approximating 

a-Quantile 







N = 


o 
1 — 1 


= 15, 


y = 


100 






TIME 


a : 


.10 


.25 


.50 


. 75 


.90 


.95 


.99 


0.01 


CLT 


0.00 


.12 


.50 


.89 


1.0 


1.0 


1.0 




HT 


0.00 


.10 


. 36 


.68 


.91 


.97 


1.0 


0.025 


CLT 


0.04 


.21 


.49 


. 81 


.97 


.99 


1.0 




HT 


0.00 


.09 


. 33 


.67 


.91 


.98 


1.0 


0.05 


CLT 


0.08 


.23 


.49 


. 77 


. 94 


.98 


1.0 




HT 


0.0 


.07 


.30 


.65 


.92 


.97 


1.0 


0.10 


CLT 


0.10 


. 24 


.48 


. 76 


.92 


.97 


1.0 




HT 


0.00 


.04 


. 24 


.60 


. 89 


.97 


1.0 



26 



TABLE 5 



Simulated Probability (Relative Frequency) that the 
Response Time is Less than or Equal to the Approximating 

a-Quantile 







N 


= 10, 


X = 25, 


U = 


100 






TIME 


a : 


. 10 


.25 


. 50 


. 75 


.90 


.95 


.99 


0.01 


CLT 


.04 


. 15 


.45 


.85 


1.0 


1.0 


1.0 




HT 


. 11 


. 22 


.43 


.72 


.91 


.98 


1.0 


0.025 


CLT 


.08 


. 19 


.45 


.80 


.98 


1.0 


1.0 




HT 


.10 


. 22 


.44 


. 73 


.92 


.98 


1.0 


0.05 


CLT 


.09 


. 22 


.45 


. 77 


.95 


.99 


1.0 




HT 


.11 


.22 


.44 


. 71 


.91 


.97 


1.0 


0.10 


CLT 


. 10 


. 24 


.47 


. 76 


.93 


.98 


1.0 




HT 


.11 


.23 


.45 


. 70 


. 89 


.95 


1.0 
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TABLE 6 



Simulated Probability (Relative Frequency) that the 



Response 


Time is 


Less 

a- 


than or 
-Quantile 


Equal 


to the 


Approximating 






N = 


25 , 


X = 10 , 


y = 100 






TIME 


a : 


.10 


.25 


.50 


.75 


.90 


.95 


.99 


0.01 


CLT 


.03 


.15 


.46 


. 86 


.99 


1.0 


1.0 




HT 


.11 


.23 


.46 


.73 


.91 


.97 


1.0 


0.025 


CLT 


.07 


. 19 


.46 


. 80 


.96 


.99 


1.0 




HT 


.11 


.23 


.46 


.73 


.91 


.97 


1.0 


0.05 


CLT 


.09 


.21 


.46 


.78 


.94 


.98 


1.0 




HT 


.11 


.24 


.46 


. 73 


.91 


.96 


1.0 


0.0625 


CLT 


.09 


.22 


.47 


. 77 


.93 


.98 


1.0 




HT 


.11 


.24 


.47 


.73 


.90 


.96 


1.0 
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5.2 Simulation Results for Markovian Model with Two-Terminal 
Types . 

In this subsection we describe the results of a simulation 
of the general K-type Markovian model of Section 4 , in the 
case of K = 2 sets of terminals. As in Section 4 let Xk(w) 
represent the number of other jobs of type i being processed when 
the tagged job has acquired exactly w units of processing. 

As before the response time for the tagged job requiring T 
units of work is 

T 

R(T) = / [X, (w) + X_(w) + 1] dw . 

0 

The process { ( (w) , (w)) ; w >_ 0 } is Markovian. Hence 
again R(T) satisfies a central limit theorem as T -* ». 

The normal approximation for the distribution of R(T) result- 
ing from the central limit theorem will again be referred to as 
CLT . 

The mean term (m^ n-m^JT for the heavy traffic approximation 
was computed by solving the system of equations (4.20) for 
m^ and • The variance term for the approximation was computed 
by solving the system of stochastic differential equations (4.23) 
as detailed in Arnold [Corollary (8.2.4) ] . The fundamental 
matrix (Arnold [p. 129]) was found by computing Laplace 
transforms of the system of defining differential equations 
and then inverting the solution. The approximating variance 
term was found by computing the variance of the solution of 
the s.d.e. As in the case of one-terminal type, it is a 
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linear combination of exponentials and constant terms. Its 
exact form is uninformative and will not be given here. 

Conditional response times, given the number of jobs of 
each type being processed at the time of arrival of the tagged 
job, were simulated. The tagged job was always taken to be 
a Type 1 job. For each initial condition, 300 applications 
were carried out. Sample moments and probabilities (relative 
frequencies) were computed for each initial condition giving 
conditional sample moments and probabilities (relative fre- 
quencies) . Unconditional sample moments and probabilities 
were computed in a similar manner to that in the one-terminal 
type simulation; see Pornsuriya [1984]. 

Values of the simulated means and standard deviations and 
their approximating values for R(T) for various cases in which 
= 5 and = 5 appear in Tables 7-10. Once again the CLT 
mean is the true steady state mean for R(T) . The means and 
standard deviations of R(T) for each T differ surprisingly 
little for the four cases. This suggests that perhaps the 
two-type terminal model can often be satisfactorily approxi- 
mated by a one-type model in which the arrival rate and service 
rates are the average arrival and service rates in the two-type 
model. Values of the simulated means and standard deviations 
and their approximating values for the approximate one-type 
model with N = 10 , A. = 25 and y = 75 appear in Table 11. The 
values for the approximate one-type model are acceptably 
close to those for the two-type model. Note that the quality 



of the HT approximation is generally quite good, even though 
the system sizes and ^ can hardly be called "large." 
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TABLE 7 





Simulated Means and Standard Deviations 
and Their Approximating Values 


for R ( T ) 




N x = 5, N 2 


= 5 , A 1 = 20 , A 2 = 30 , = 


: 50, p 2 = 


TIME T 




Mean 


Std. Dev. 


0.01 


Simulation 


0 .0713 
(0.0001) * 


0.0135 




CLT 


0 .0707 


0.0204 




HT 


0 .0710 


0.0132 


0.025 


Simulation 


0.1783 

(0.0003) 


0.0263 




CLT 


0.1766 


0.0322 




HT 


0 .1775 


0.0253 


0.0375 


Simulation 


0.2664 

(0.0005) 


0 .0353 




CLT 


0 . 2650 


0.0395 




HT 


0.2663 


0 .0325 


0.0 50 


Simulat ion 


0 . 3570 
(0 .0005) 


0.0414 




CLT 


0 . 3533 


0.0456 




HT 


0 . 3550 


0 .0384 


0.0625 


Simulation 


0.4439 

(0.0006) 


0.0478 




CLT 


0.4416 


0.0510 




HT 


0.4438 


0.0435 



* 

Standard error of the estimate of the mean 



100 
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TABLE 8 





Simulated Means and Standard Deviations 
and Their Approximating Values 


for R ( T) 




N 1 = 5 ' n 2 = 5 ' 


A 1 


= 30 , X 2 = 20 , ] a ± = 


100, y 2 = 


TIME T 






Mean 


Std. Dev 


0 .01 


Simulation 




0.0710 

(0.0001)* 


0.0137 




CLT 




0.0716 


0.0204 




HT 




0.0710 


0.0132 


0 .025 


Simulation 




0.1777 
(0 . 0003) 


0.0268 




CLT 




0.1789 


0 .0322 




HT 




0 .1775 


0.0253 


0.0375 


Simulat ion 




0.2672 

(0.0004) 


0.0351 




CLT 




0.2684 


0.0395 




HT 




0.2663 


0 .0325 


0.050 


Simulation 




0 . 3567 
(0.0005) 


0 .0419 




CLT 




0 . 3578 


0.0456 




HT 




0 . 3550 


0.0384 


0.0625 


Simulation 




0.4461 

(0.0006) 


0.0476 




CLT 




0.4473 


0.0510 




HT 




0 .4438 


0.0435 



Standard error of the estimate of the mean 
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TABLE 9 



Simulated Means and Standard Deviations for R(T) 
and Their Approximating Values 





N 1 = 5/ N 2 = 5, 


o 

II 

CM 

o 
1 — 1 

II 

1 — 1 


CM 

n 

ID 

CM 

II 

i — ! 


TIME T 




Mean 


Std. Dev 


0 .01 


Simulation 


0.0717 
(0.0001) * 


0 .0134 




CLT 


0.0717 


0.0237 




HT 


0.0720 


0.0133 


0.025 


Simulation 


0.1788 

(0.0004) 


0.0280 




CLT 


0 .1793 


0 .0375 




HT 


0.1799 


0 .0271 


0 .0375 


Simulation 


0.2690 

(0.0005) 


0.037 6 




CLT 


0 .2684 


0.0459 




HT 


0.2699 


0.0359 


0.0500 


Simulation 


0 . 3588 
(0.0006) 


0.0456 




CLT 


0 . 3585 


0.0530 




HT 


0 . 3599 


0.0433 


0.0625 


Simulation 


0 .4481 
(0.0007) 


0.0527 




CLT 


0.4481 


0.0529 




HT 


0 .4498 


0.0496 



* 

Standard error of the estimate of the mean 



TABLE 10 



Simulated Means and Standard Deviations for R(T) 
and Their Approximating Values 





N 1 = 5 ' N 2 = 5 ' 


o 
1 — 1 

II 

04 

-< 

O 

II 

r — 1 

r< 


P 1 = 125, y 2 


TIME T 




Mean 


Std. Dev. 


0.01 


Simulation 


0.0725 
(0.0001) * 


0.0138 




CLT 


0.0724 


0.0249 




HT 


0.0720 


0.0133 


0.025 


Simulation 


0.1816 

(0.0004) 


0.0290 




CLT 


0 . 1810 


0 .039 4 




HT 


0.1799 


0.0271 


0.0375 


Simulation 


0.2717 

(0.0005) 


0 .0389 




CLT 


0.2714 


0.0482 




HT 


0 .2699 


0.0359 


0.0500 


Simulation 


0 . 3622 
(0.0006) 


0.0472 




CLT 


0.3619 


0.0557 




HT 


0.3599 


0.0433 


0.0625 


Simulation 


0.4522 

(0.0008) 


0.0557 




CLT 


0.4524 


0.0623 




HT 


0.4498 


0 . 0496 



*Standard error of the estimate of the mean 
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TABLE 11 





Simulated Means and Standard Deviations 
and Their Approximating Values 
for the One-Type Model 


for R ( T) 




N = 


10, 


A = 25, y = 75 




TIME T 






Mean 


Std. Dev. 


0 . 1 


Simulation 




0 .0701 
(0.0002) * 


0.0139 




CLT 




0 .0701 


0.0206 




HT 




0.0700 


0.0135 


0.025 


Simulation 




0.1766 

(0.0005) 


0.0263 




CLT 




0 . 1752 


0.0325 




HT 




0 . 1750 


0.0258 


0.0375 


Simulation 




0.2620 

(0.0008) 


0.0361 




CLT 




0 . 2628 


0 .039 8 




HT 




0 .2625 


0.0330 


0 .0500 


Simulation 




0 . 3501 
(0.0009) 


0.0430 




CLT 




0 . 3504 


0.0460 




HT 




0 .3500 


0.0390 


0 .0625 


S imulation 




0 .4388 
(0.0011) 


0.0478 




CLT 




0.4380 


0.0514 




HT 




0.4375 


0.0441 



★ 

Standard error of the estimate of the mean 
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To assess the quality of the normal approximation to the 
distribution of R(T) for the two- type model, the a-quantiles 
for each two-type approximating normal distribution were 
computed. The relative frequency of being less than or equal 
to each a-quantile was then computed, using the simulated 
data. The results appear in Tables 12-15. From the heavy 
traffic approximation to the mean it follows that approximately 
7 jobs are being processed with the tagged job. Thus, all 
the cases considered are really moderate traffic cases. The 
tables indicated that the HT approximation tends to describe 
the quantiles better than does the CLT. However, as is 
expected, the CLT improves with larger T. 
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TABLE 12 



Simulated Probability (Relative Frequency) that the 
Response Time is Less than or Equal to the Approximating 

a-Quantiles 





LO 

II 

f — 1 


N 2 = 5 


II 

i — 1 


20, X 2 


II 

u> 

o 

TT 


x = 50, 


U 2 


100 


TIME T 


a : 


0.10 


0.25 


0.50 


0 . 75 


0.90 


0 .95 


0.99 


0 .010 


CLT 


0.041 


0.137 


0.433 


0 . 838 


0.998 


1.0 


1.0 




HT 


0.109 


0.221 


0.443 


0.719 


0.915 


0.977 


1.0 


0.0250 


CLT 


0 .062 


0.170 


0.425 


0 . 765 


0.969 


0.997 


1.0 




HT 


0.103 


0.223 


0.438 


0 . 708 


0.911 


0.975 


1 . 0 


0.0375 


CLT 


0 .079 


0.191 


0 .429 


0 . 754 


0.952 


0.991 


1.0 




HT 


0 . 118 


0.229 


0.447 


0 . 713 


0 .9 05 


0.969 


0.999 


0.0500 


CLT 


0 .075 


0.181 


0.417 


0 . 724 


0 . 939 


0.986 


1.0 




HT 


0.106 


0 . 222 


0.431 


0 . 693 


0 . 89 7 


0.965 


0.998 


0.0625 


CLT 


0.081 


0.198 


0.436 


0.745 


0.932 


0.980 


1.0 




HT 


0.115 


0 .238 


0.451 


0.719 


0 . 898 


0 .960 


0.997 



TABLE 13 



Simulated Probability (Relative Frequency) that the 
Response Time is Less than or Equal to the Approximating 

a-Quantiles 





N 1 = 


II 

CN 

ID 


II 

i — 1 

ID 


30, A 2 


o 

CN 

II 


v 1 = ioo 


y 2 = 


50 


TIME T 


a : 


0.10 


0 .25 


0.50 


0.75 


0.90 


0.95 


0.99 


0.0100 


CLT 


0.049 


0.159 


0.457 


0 . 869 


0.998 


1 . 0 


1.0 




HT 


0.114 


0 . 229 


0.442 


0.718 


0.919 


0 .979 


1.0 


0.0250 


CLT 


0.078 


0 .194 


0 .462 


0.809 


0.978 


0.997 


1.0 




HT 


0 .113 


0 .228 


0.442 


0.717 


0.916 


0.973 


0.999 


0.0375 


CLT 


0.088 


0 .204 


0.455 


0 . 787 


0.962 


0.995 


1.0 




HT 


0.111 


0 .224 


0.432 


0 . 702 


0.909 


0.965 


0.999 


0.0500 


CLT 


0.088 


0.212 


0 . 459 


0 .771 


0 .955 


0 .9 89 


1.0 




HT 


0.111 


0.226 


0.4 34 


0.696 


0 . 89 8 


0.961 


0.998 


0.0625 


CLT 


0.093 


0 .212 


0.465 


0 . 773 


0.945 


0.986 


1.0 




HT 


0 .110 


0.222 


0 .432 


0 . 705 


0 . 889 


0 .954 


0.997 
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TABLE 14 



Simulated Probability (Relative Frequency) that the 
Response Time is Less than or Equal to the Approximating 

a-Quantiles 





II 

i — 1 


5, N 2 


* — 1 

m 

II 


= 10, 


= 40 

2 


' y l “ 


25, y 2 


= 125 


TIME T 


a 


0.10 


0 . 25 


0.50 


0.75 


0.90 


0.95 


0.99 


0 . 010 


CLT 


0.027 


0.119 


0.447 


0.906 


1.0 


1.0 


1.0 




HT 


0 .110 


0.230 


0.458 


0 . 737 


0 .9 31 


0.982 


1.0 


0.0250 


CLT 


0.061 


0.178 


0 .451 


0 . 824 


0.991 


0.999 


1.0 




HT 


0.121 


0.236 


0.464 


0 . 734 


0.930 


0.985 


1.0 


0.0375 


CLT 


0.068 


0.191 


0 . 441 


0 . 790 


0.979 


0.999 


1.0 




HT 


0 . 117 


0.239 


0 .451 


0.725 


0.927 


0.983 


1.0 


0.0500 


CLT 


0.081 


0 .189 


0.430 


0 . 776 


0.975 


0.997 


1.0 




HT 


0.117 


0.232 


0.444 


0 . 729 


0.935 


0.985 


1.0 


0.0625 


CLT 


0.085 


0.199 


0.442 


0 .772 


0.960 


0.995 


1.0 




HT 


0 .121 


0 . 237 


0.456 


0.737 


0.923 


0 .979 


1.0 
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TABLE 15 



Simulated Probability (Relative Frequency) that the 
Response Time is Less than or Equal to the Approximating 

a-Quantiles 





N 1 = 


5, N 2 


i — 1 

ID 

II 


= 40, 


o 
• — 1 

II 

CM 


' P 1 " 


125, u 2 


= 25 


TIME T 


a : 


0 . 10 


0.25 


0 . 50 


0 . 75 


0.90 


0.95 


0.99 


0.0100 


CLT 


0 .024 


0.116 


0.445 


0.907 


1.0 


1.0 


1.0 




HT 


0.110 


0 . 228 


0.433 


0 . 710 


0.905 


0.9 73 


1.0 


0.0250 


CLT 


0.054 


0 .164 


0.435 


0 . 820 


0.993 


1.0 


1.0 




HT 


0.112 


0 . 219 


0.422 


0 .682 


0.901 


0.967 


1.0 


0.0375 


CLT 


0.070 


0.186 


0.428 


0.79 5 


0.981 


0.999 


1 . 0 




HT 


0.113 


0.224 


0.414 


0 . 693 


0. 89 9 


0.967 


0.999 


0 . 050 


CLT 


0.075 


0.182 


0.439 


0 . 784 


0.972 


0.998 


1.0 




HT 


0.110 


0.216 


0.422 


0 .689 


0.89 6 


0.969 


1.0 


0.0625 


CLT 


0.083 


0.200 


0.443 


0 . 765 


0.959 


0.997 


1.0 




HT 


0 . 118 


0.222 


0.427 


0.683 


0 . 891 


0.958 


1.0 
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APPENDIX 



HEAVY TRAFFIC APPROXIMATION BY CONVERGENCE- 
0 F- S EMI GRO UP S METHODOLOGY 

The purpose of this appendix is to outline a ma themazical 
framework upon which the heavy traffic approximations of this 
paper may be rigorously based. The approach is to use an 
analytical theory of convergence of semigroups of operators 
apparently first, applied to queueing problems by Burman (19 79) 
m a regrettably unpublished thesis. See also Lehoczky and 
Gaver (1981) where the technique is used to obtain results 
concerning a data-voice traffic sharing multichannel system. 

The theory of semigroups of operators is introduced in Feller 
19~1), and detailed in Dynkin (1965); the convergence ideas 
are discussed in Trotter (1974) and Kato (1976) . The basic 
notion is that the state variable of a process, say the work 
time process of Section 4.1, {x^(w;N), w >_ 0 } , is one of a 
sequence of birth and death Markov processes indexed by system 
size N. Given such a sequence of Markov' processes, <f X^(w;N)}> , 
each with its appropriate state space, S^ , it is desired to 
show that the corresponding sequence of probability transition 
functions converges to that of some limiting process that has 
state SDace S ; in the present case S = TR . The limiting 

-co OO -f- 

process under the normalization of {X(w;N)} chosen will be a 
particular diffusion process, namely, in the present heavy 
traffic situation the multivariate Ornstein-Uhlenbeck . 

The Trctter-Kato theory of convergence deals with the 
convergence of infinitesimal operators A. of the normalized 
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processes {X^(w) ;N}: if A^f -*• A^f in sup norm for a suitable 

class of test functions (e.g. f(z): (R ->■ H m-times continuous ly 

differentiable, m >_ 3 , that vanish identically outside a bounded 

If 

subset of TR and further such that the functions, f, together with their 
first and second derivatives do not increase faster than some 
fixed power of z) it can be concluded that the semigroups 
converge, and hence the Markov probability transition functions 
themselves converge. 

We now proceed with the formal calculation of the limiting 
generator for our normalized process. Invoke (4.4) so 

.. X N (w) -N.m.(w) X N (w) - NJl.m.(w) 

Z N (w) = = -i ^ . (A- 1 ) 

1 /N. /ITT /N 

l l 



By definition, for z = ( z, , z_ , . . . , z T .) and f in the above class 

1 2 K 



A M f(£) = lim{E [f ( Z_ N (w+A) ) | Z_ N (w) = z] - f(£)}j 
A^O 



(a- 2 ; 



Given Z^(w) = z^, and C^(w,w+A) represents the change in X^(w) , 

M X^(w) + C. (w,w+A) - N.m. (w+A) 

Z W (w+A ) = — 

1 /N~ 

l 



C, (w) 



l 



+ z. - /NT m!(w)A+ o(A) . 

i ii 



(A- 3 ) 



Consequently, for z such that z. > -/N7m. (w) A - — — 

1 11 /N. 

l 
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lim[ I {f(z,,...,z. + — /r~/Nm ! (w) A , . . . , z T J A 

A-0 i=l 1 1 /I“/N 11 K 

l 



+ f ( z 1 



z. - — - Si . /Nm ! (w) A , . . . , z ) y . ( N ) A } 

1 /TVn 11 K 1 

1 



+ f ( z - /T7/Nm ' (w) A , . . . , z (w ) A ) 

-L J- -L K K K 

K 

X (1 - £ { (A. (N) + y . (N) ) A}) 

i = l 1 1 

* r z ^ f * * * * ^ ^?) o ( A ) ] ~ 



where for simplicity (and generality) we abbreviate 



_ K 

A-(N) = A • [N . - n . rn . (w) -/N.z.J ( j (N .m . (w) +/N . z . ) ) 

1 1111 i i J J 3 3 



Si . K SIT. 

( A . N) [ l . ( 1-m . (w) ) - — z.](N T (£.m.(w) + 1 z . 

11 SN 1 j = l 33 SW 3 



and 



Sl~ 

y . (N) = y . (N.m (w) + /Nz.) ~ y . N ( £ • m . (w) + — - z.) . 

l ill l ill /n 1 



Upon passage to the limit via Taylor series expansion it is 



seen that 



i (N) A 



(A-4) 



( A-5 ) 

) ) 

(A-6) 
(A- 7 ) 
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1 



K 



l {f (z, , . 

i=l 1 



. Z . + 

1 



y/1 . /N 

i 



. , Z ) A . (N) 
K l 



+ f ( , . . . , - 



, . . . , z ) y . ( N ) } 

/£ . /N K 



K 



+ f ( z 1 , . . . , . 



'V [ ~ 



1 



I ( A . (N) + y . (N) ) ] 

= 1 3 3 



K 

i f; 

j =i j 



'T~ /N m ’ (w) 

j j 



(A- 8 ) 



Note that no specific normalization has been utilized up to 
this point. Now, however, invoke the HTN of (4.6) and utilize 
(A-6) and (A-7) specifically; allowing N to become large, 



K A (N) - y . (N) , 

V (z) ' /r/s 

l 



A, (N) + u i (N) 



2. • N 

l 



- f / 1 . /N m ! (w) } 
z . 1 1 

l 



( A- 9 ) 



K 



K 



(w) - y^/if~iTK(w) 



y f ’ (/N[a: / 2 . ( 1 -m. (w) ) y l.m. 

1=1 Z i 111 j=l 3 3 

K K 

- /£ . m ! (w) ] + A ! /T7 ( 1-m . (w )) ( T /2.z.)-A!z. T 2.m.(w) 
11 11 1 ^ = i 3 3 11 - = i 1 1 



- y.z. + o ( N ^ //2 ) } + ^ [ f" { A .' ( 1-m . (w ) ) 

ii 2 . c . z • i i 



K 

l 

i=l “i 



K 



2 m (w) + y . m . 

j=l j D 11 



(w) + 0(N~ 1/2 ) } 
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Choose the functions nr so that they satisfy the system of 
differential equations (4.5) . Let N ->• 00 . Then for f in the 
above class of functions, the operator converges to yield 

K K K 

A co f(2) = I f z UJ/I-U-m.Cw)) l /T7z. - A ! z . £ £ m. (w)-y.z.} 

i=l z i 1 1 1 j=l : J 1 1 j = l ^ 3 11 

, K K 

+ y l f" { A ! (1-m . (w) ) ( l £.m.(w)) + y.m.(w)} 

Z i=l z i 1 1 j=l 3 J 11 

The operator A^ is the infinitesimal operator of the diffusion 
whose stochastic differential equation is (4.6) [cf. Arnold 
(1974), page 152]. The Trotter-Kato theorem can now be applied 
to assert that the semigroups converge [cf. Burman (1979)]. 
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